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Characterization of periodic cavitation in an optical tweezer 
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Microscopic vapor explosions or cavitation bubbles can be generated periodically in an optical 
tweezer with a microparticle that partially absorbs at the trapping laser wavelength. In this work 
we measure the size distribution and the production rate of cavitation bubbles for microparticles 
with a diameter of 3 /im using high speed video recording and a fast photodiode. We find that there 
is a lower bound for the maximum bubble radius Rmax ~ 2 /i m which can be explained in terms of 
the microparticle size. More than 94% of the measured Rmax are in the range between 2 and 6 /im, 
while the same percentage of the measured individual frequencies fi or production rates are between 
10 and 200 Hz. The photodiode signal yields an upper bound for the lifetime of the bubbles, which 
is at most twice the value predicted by the Rayleigh equation. We also report empirical relations 
between Rmax , fi and the bubble lifetimes. 


I. INTRODUCTION 

Laser-induced cavitation involves the fast deposition 
of laser energy in a liquid, resulting in an explosion that 
creates a transient vapor bubble that expands reaching 
a maximum radius Rmax and then collapses under the 
static pressure of the surrounding liquid [1]. These bub¬ 
bles have been used in many different applications where 
fast actuation of liquid or where impulsive forces are re¬ 
quired at a microscopic length scale. For example, in high 
speed microfluidics [2] and micro-pumps [3, cell mem¬ 
brane permeabilization Si. and cell lysis [fi, red blood 
cell stretching and poration [7, 8], malaria detection [9|, 
bending of carbon nanotubes and nanowires mm- 

In the last few years there has been a lot of inter¬ 
est in the production of cavitation bubbles through pho- 
tothermal processes in nanomaterials [3, M, Il2l Il4j result¬ 
ing in bubbles with sizes on the order of several hun¬ 
dreds of nanometers. These transient nanobubbles have 
shown great potential for medical treatments El- op- 
tical tweezers or highly focused continuous (CW) lasers 
have been used to create micron sized cavitation bub¬ 
bles in liquids that are transparent to the laser wave¬ 
length pi| [13 • Dholakia and coworkers [l6| trapped a 
polystyrene nanoparticle that was later irradiated with 
a nanosecond laser pulse creating optical breakdown at 
the trapped nanoparticle. Another approach used cus¬ 
tom made microparticles filled with absorbing dye that 
was rapidly vaporized creating multiple cycles of bubble 
expansion and collapse in a timescale of a few microsec¬ 
onds Q3. 

Recently we showed that a microparticle that partially 
absorbs at the wavelength of the trapping beam in an 
optical tweezer can periodically generate microscopic ex¬ 
plosions [HI- The particle is simultaneously heated and 
pulled towards the waist of the trapping beam where it 
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superheats the liquid and creates a transient vapor ex¬ 
plosion that pushes the particle below the waist, where 
the cycle restarts (Fig. la). However, in that study we 
could not establish the bubble size distribution of Rmax 
and production rate, which are important for potential 
applications. In the present work we measure hundreds 
of bubbles with two different methods: fast video record¬ 
ing and a fast photodiode which dramatically increases 
the time resolution. We find that there is a minimum 
bubble size of Rmax ~ 2/im and that most of the bub¬ 
bles have an R m ax between 2 and 6 /im which is a fairly 
narrow range. 

II. EXPERIMENT 

The experimental setup is a near IR optical tweezer 
with a wavelength of 975 nm focused by a 100X/1.25 NA 
microscope objective. The microparticles are standard 
magnetic beads (Bangs labs) with a mean diameter of 
3.16 /im immersed in water. In the experiments, the mi¬ 
croparticles interact individually with the focused laser. 
The liquid sample is placed between two No. 1 micro¬ 
scope glass cover slips (thickness 0.13 — 0.16 mm) sepa¬ 
rated with spacers with a height ^100 /im. The sam¬ 
ple is then sealed with nail polish. Some experiments 
were done without sealing the sample and the results 
were consistent with those with the sealed chamber. A 
schematic of the experimental setup is depicted in Fig. 
lb. The transmitted laser power by the microscope ob¬ 
jective is between 50 and 62 mW. The trapping laser light 
is collected with the microscope condenser and a lens into 
a fast photodiode (175ps rise, Alphalas) which replaces 
the quadrant detector in a typical optical tweezer setup. 
The photodiode is connected to an oscilloscope (1GHz, 
LeCroy 610Zi). The sample is illuminated by an LED 
and imaged into a high speed video recorder (Photron 
SA 1.1). 

The oscilloscope and the high speed video recorder 
are triggered with a delay generator (SRS DG535). The 
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FIG. 1: (a) Periodic cavitation in an optical tweezer, (b) Experimental setup, (c) Selected frames corresponding to (a) from a 
high speed recording at 300,000 fps. (d) Light scattered by the microbead and the bubble detected with the photodiode (black 
line). Particle position p(t) in the XY plane (red line), p(t) = y/(X(t) — Xq) 2 + (Y(t) — Yo) 2 where Xq, Yq are the coordinates 
of the beam center. 


high speed video recorder captures between 100,000 and 
150,000 frames at 300,000 frames per second (fps). The 
oscilloscope records 25 million points in a time span of 
one second which yields a time step At of 40 ns, sufficient 
to measure the width of the signal of the scattered light 
during cavitation. Typical images extracted from a high 
speed video recording of a single cycle are shown in Fig. 
lc. First the particle is below the focus of the microscope 
objective which is why it is blurred, then it is pulled to¬ 
wards the geometrical focus where the particle appears 
sharper. The cavitation bubble appears in one frame and 
in the next (3.3 /is later) the bubble has vanished and the 
particle is below the focus of the objective. 

A sample of an oscilloscope trace for the scattered light 
collected by the photodiode is shown in Fig. Id (black 
line, 0-0.5 s). The bubbles appear at the sharp transi¬ 
tions, labeled as the i-th bubble which is followed by a 
time interval Ti (cycle period) that depends on the size 
of the bubble. The frequency of that cycle is fi = 1/T*. 
The red line corresponds to the distance of the micropar¬ 
ticle in the XY plane to the center of the trap (located 
at (X 0 , Yo)) where the laser beam is focused. 


III. RESULTS 


Figure 2a (Left) shows the zoomed photodiode signal 
for selected cavitation events and the corresponding im¬ 
ages for the bubble and particle extracted from the high 
speed videos. We observe that there is a change in the 
voltage level which corresponds to the particle being close 
to the waist of the beam (before cavitation) and after the 
particle has been pushed by the bubble below the waist. 
The transition signal in between captures the large parti¬ 
cle displacement and the bubble dynamics which is why 
it can have different shapes (peak or a valley) depend¬ 
ing on the spatial location where cavitation occurs and 
the direction in which the particle is pushed due to rapid 
changes in Mie scattering by both particle and bubble. 


Usually the transition starts with a small peak followed 
by a larger one, hence it is not straighforward to locate 
the part that corresponds to the bubble. The width T p d 
of the transition is defined by numerically differentiating 
(finite difference) the signal, taking the absolute value 
and measuring the time between the two largest peaks 
(insets in Fig. 2a Left). One more condition is that the 
two maximums (in the derivative) have to span the whole 
transition (including the smaller first peak), in order to 
get the maximum width and hence an upper bound to 
the bubble lifetime T p d measured with the photodiode. 
This criteria yields 679 events. 

We measure the maximum bubble radius R m ax from a 
single frame in the high speed videos where the bubble 
appears blurred (Fig. 2a Right). The error in the mea¬ 
surement of Rmax is ±0.22 pm estimated by varying the 
gray level thresholds that define the edge of the bubble. 
We were able to measure 179 bubbles where the image 
had good signal to noise ratio. To estimate the bubble 
lifetime (the time it takes the bubble to fully expand and 
collapse) we assume spherical symmetric bubble dynam¬ 
ics [UJ, [20j . The spherical bubble lifetime is T osc which is 
twice the collapse time Tc that can be calculated using 
the Rayleigh formula 0: 


1.82 Rr, 


Pi 


PO -Pv 


(1) 


where pi is the density of the liquid (water), p Q the at¬ 
mospheric pressure (0.1 MPa) and p v the vapor pres¬ 
sure at room temperature (2330 Pa at 20°C ). However 
there is a loss in the accuracy of the Rayleigh formula 
for nanobubbles 0 (since it does not take into account 
surface tension). In this study the maximum deviation 
should be at most a factor of 1.2 0 for the smallest bub¬ 
bles Rmax <2.7 pm. Hence we use T osc to get a lower 
bound for the bubble lifetime. The spherical shape as¬ 
sumption is a strong one but it has been used successfully 
to model microparticles accelerated by cavitation bub¬ 
bles created at the microparticle surface [22|. Further- 
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FIG. 2: (a) Left column: Photodiode signals during cavitation events for different bubble sizes. Inset: Absolute value of the 
numerical first derivative in the same temporal range, (a) Right column: Frame where the bubble corresponding to the event 
in the left appears, (b) Data for calculated spherical oscillation time (from measured Rmax) as a function of the width of the 
photodiode signal (116 events). 


more, some experimental measurements on the dynamics 
of photothermal bubbles show symmetric dynamics [5]. 

For some events (116) we were able to measure si¬ 
multaneously T p d and Rmax • The plot of the calcu¬ 
lated oscillation time T osc (extracted from the measured 
Rmax) as a function of the width of the photodiode sig¬ 
nal T p d is in Fig. 2b. The bubble lifetime predicted 
by the Rayleigh formula is proportional to the measured 
width of the photodiode signal during the transition: 
Tosc = 0.45 Tpd + 121.3 ns. The width of T p d is an upper 
bound to the lifetime, so Fig. 2b shows that the lifetime 
is at most two times that predicted by Rayleigh equa¬ 
tion for T osc [Eq. (1)]. Deviations could come from non 
spherical dynamics, compressibility and other photother¬ 
mal processes that have been neglected. 

Figure 3(a-c) depicts the data extracted from the high 
speed videos (179 events). The histograms in in Figs. 3a 
and 3b correspond to the measured Rmax and the indi¬ 
vidual frequencies fi respectively. The largest bars are: 
31 events where R m ax has values between 2.95 — 3.22 /i m 
and 38 events where fi has a value between 24.1—35.7 Hz. 
There are 13 events with R ma x between 1.85-2.12 fi m, 
one event with an R m ax of 1.6 /im and no events with 
a smaller R m ax • The individual cycle frequencies fi are 
plotted in Fig. 3c as a function of the measured Rmax • 
The red solid line is a fit to an exponential decaying 
function fi = a + 6exp (—ci? macc ), where a = 14.7 Hz, 
b = 1434.8 Hz and c = 1.0412 /im -1 . The fit is rea¬ 
sonable up to frequencies of about 200 Hz. This could 
be useful to estimate bubble size by just measuring the 
frequency. 

The photodiode data (679 events) are shown in Figs. 
3d and 3e which are histograms for T p d and fi respec¬ 
tively. The largest bars in the histograms are 106 events 
with Tpd between 1010-1090 ns and 106 with fi between 
63.9-75.5 Hz. The data in Fig. 3d can be used to get 
Rmax by first using the linear relation between T osc and 
Tpd and then solving for R m ax in [Eq. (1)]. There are 106 
events with R m ax between 3.15 — 3.34 /im (Fig. 3d). Sur¬ 
prisingly we get the same lower bound for Rmax ~ 2 /im 


as in Fig. 3a, considering that the time resolution with 
the photodiode is much greater than with the video cam¬ 
era and can resolve very short T p d- There is 1 event 
corresponding to 1.85 /im, 3 events at 1.95 — 2.15 /im 
and 17 at 2.15 — 2.35 /im. 


IV. DISCUSSION AND ANALYSIS 

The parameters that describe the interaction between 
a spherical microparticle and a cavitation bubble are the 
microbead radius R p , the maximum bubble radius Rmax , 
the initial separation between the microparticle and the 
origin of the bubble S' and the displacement of the par¬ 
ticle Ax' 0®. The interaction is described by a mas- 
tercurve for Ax as a function of Sx where S = S'/Rmax 
and Ax = Ax' /(2 Rp). After the bubble collapses the mi¬ 
croparticle can stop farther away from the starting point 
(repulsion Ax > 0 for S < 0.7), closer to the starting 
position (attraction Ax < 0 for 0.7 < S < 3.7) or at the 
same distance (neutral Ax = 0 for S = 0.7). 

During the bubble expansion the particle is displaced 
in a timescale of microseconds (repulsion Ax > 0, S < 
0.7), while the return to the beam waist with optical 
forces is much slower in the tens of milliseconds. Hence 
it is expected that the bubbles with the smaller R ma x 
will induce a smaller displacement and a shorter time 
back to the position of the next vapor explosion. Since 
the bubble starts at the surface of the microparticle then 
S' = R p = 1.58 /im. The equilibrium position for S = 0.7 
would correspond to the maximum frequency since es¬ 
sentially the particle would end at the same spot after 
the cavitation bubble collapses where it should be su¬ 
perheated again in a very short time producing another 
bubble. For S = 0.7, Rmax = 2.26 fi m, which is con¬ 
sistent with the smallest measured Rmax (Figs. 3a, 3d) 
which correspond to the largest frequencies (Fig. 3c). 

The width in the histograms for fi presented in Figs. 
3b and 3e is larger than that for R m ax and T p d . This can 
be explained considering the random direction in which 
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FIG. 3: (a) Histogram for measured Rmax (video data 179 events), (b) Histogram for measured /* (video data), (c) Cycle 

frequency fi as a function of maximum bubble radius Rmax (video data), (d) Histogram for measured T p d (photodiode data 
679 events), (e) Histogram for measured fi (photodiode data). 


the particle is pushed which results in different optical 
forces (oc V/, with / the intensity) and hence different 
speeds. 

The distribution of Rmax in Figs. 3a and 3d is consis¬ 
tent with the previous measurements [l8| where only 40 
bubbles were measured and smaller bubbles were not de¬ 
tected due to insufficient signal to noise caused by poorer 
illumination of the sample. Here the measured frequen¬ 
cies lie in a larger interval than that reported in |18| (23 
to 48 Hz) because in that study the frequencies were ex¬ 
tracted from the largest peaks on the Fourier spectrum 
of the microparticle dynamics. 

V. CONCLUSION 

We have measured distribution of Rmax which is sur¬ 
prisingly narrow with most values between 2-6 /im. The 
size is comparable with bubbles created with focused fem¬ 
tosecond laser pulses uu and some photothermal bub¬ 
bles in nanomaterials J5[. The measured cutoff or lower 
bound for the bubble sizes of ^ 2 fi m a gree s with the 
neutral displacement regime reported in [23|. We also 


found that the production rate of the cavitation bubbles 
depends in the bubble size as expected, since a larger 
bubble results in a larger displacement for the micropar¬ 
ticle and a slower return back to the waist of the trapping 
beam. The width of the photodiode signal yielded upper 
bounds to the lifetime of the bubbles which are at most 
twice those predicted by the Rayleigh formula for T osc . 
The empirical relations between R m ax , T osc , T p d and fi 
could be useful to estimate Rmax and the lifetime when a 
fast video recorder is not available. Further studies in the 
shape and dynamics of the bubbles could include ultra- 
high speed imaging with temporal resolution of tens of 
nanoseconds in order to measure the full dynamics and 
shape of the bubbles. 
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